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1.0 INTRODUCTION 


Graphics is one of the important tools in engineering analysis and 
design. However, plotting routines that generate output on high 
quality laser printers normally come in graphics packages, which 
tend to be expensive and system dependent. These factors become 
important for small computer systems or desktop computers, 
especially when only some form of a simple plotting routine is 
sufficient. With the Postscript language becoming popular, there are 
more and more Postscript laser printers now available. Simple, 
versatile, low cost plotting routines that can generate output on 
high quality laser printers, are needed and standard Fortran language 
plotting routines using output in Postscript language seems logical. 
The purpose of this report is to explain the use of two simple 
Fortran plotting routines that generate output in Postscript 
language. 
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2.0 OVERVIEW OF THE PROGRAMS 


The listing of two Fortran plotting programs, PSPLOTXY.FOR and 
PSTWINPLOT.FOR, is in the appendices. The first program generates 
plots in X-Y format. A sample output is show in figure 1. The 
second program generates two plots on the same page (called 
twinplot). A possible and popular application for the second program 
is to display both phase and magnitude for a set of complete data on 
the same page. A sample plot in this format is shown in figure 2. 

These two routines can generate multiple curves (a maximum of 
4) on the same plot with different line styles. Users have options to 
label curves, label the axes, and title the plot In addition, two lines 
of comments can be included. Furthermore, the routine can 
automatically list the plot date and the time the plots are 
generated. However, this option is system dependent and it makes 
the programs nonportable. Users have the option to include this 
feature at the cost of code portability. This is described in section 
5.0 in detail. 

The two programs and their inputs are similar, thus the 
instructions for .their use can be combined. 
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3.0 PROGRAM STRUCTURES 


The main program is very simple. It makes two calls to the 
subroutines. The first call is to subroutine DATAIN to read data to 
be plotted. The second call is to subroutine PPLOT to generate the 
plot. The calls to these subroutines are repeated for multiple plots 

The subroutine that reads in the data requires that the data be 
formatted in one of two specific formats. These formats are 
explained in section 4.0. 

The second subroutine generates all the plots, labels, etc. The 
explanation of this subroutine requires some knowledge of the 
Postscript language [1], which is beyond the scope of this report, and 
will be skipped. 
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4.0 INPUT FILE FORMAT 


The subroutine that reads data to be plotted is named DATAIN. 
The data should be stored in a file in one of the two following 
formats: 

Format 1 : 

NPTS, XMIN, DX 
AM(1), PH(1) 

AM(2), PH(2) 

AM(3), PH(3) 


The first line represents the number of data points, NPTS, the 
starting X value, XMIN, and the increment DX between successive 
values of X. This format assumes data are equaliy spaced in the X 
coordinate with XMIN being X(1). 

The input lines following the first represent the amplitude, 
AM, and phase, PH, values of the data. AM and PH are values in the 
vertical coordinates. AM and PH represent Amplitude (and Phase 
arrays, respectively (although the data may be something else if the 
users wish). For XY plots (PSPLOTXY FOR), only the AM array is used 
in plotting. However, the PH array is still required in the input to 
maintain input file compatibility for both programs PSPLOTXY. FOR 
and PSTWINPLOT.FOR. 

Form at 2: 

NPTS 

X(1), AM(1), PH(1) 

X(2), AM(2), PH(2) 

X(3), AM(3), PH(3) 


i 
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The first line represents ^hfe number of points, N.PTS. 

The following lines represent the X values for the horizontal axis 
and the AM and PH values for the vertical axis. For XY plots 
(PSPLOYXY FOR), only the X and AM arrays are used. However, the 
users must specify dummy PH values for the data to be read 
properly. 
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5.0 INTERACTIVE INPUT 


Program users compile and run codes just like any other Fortran 
program. The programs prompt the users for input. The procedure is 
very self-explanatory. For completeness, the inputs are explained 
below: 

Input 1 : NPLOTS 


NPLOTS is the number of curves to be plotted on the same plot; the 
maximum number of curves to be plotted is four. 

Input 2 : IDEF 

IDEF is the axis label indicator. 

I 

IDEF = 0 For automatic axis labeling of Magnitude vs. 

Frequency, 

= 1 For automatic axis labeling of Magnitude vs Angle 
plots, 

* 2 For the program to prompt users to input axis labels, 

= 9 To change plot size and prompt for axis labels. 

The default plot size is 8.0 X 5.5 inches for XY plot (figure 1) and 
5.688 x 3.438 inches for each plot on the twinplot (figure 2). 

Input 3 : I FI 

IFI indicates which format the data is in, as explained previously 
Input 4 : INFILE 

Input data file name. The file must have the format specified by 
Input 3. 
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Input 4a : 


If IDEF = 2 users are prompted for axis labels for horizontal and 
vertical axis. 

If IDEF = 9 users are prompted for a new plot size and axis labels. 
Input 5 : XMIN, XMAX, XINC 

After displaying the minimum and maximum values of X array, users 
are prompted for the lower (XMIN) and upper (XMAX) limits of plot, 
and the X increment (XINC) between adjacent tic marks for the 
horizontal axis. 

Input 6 : YMIN, YMAX, YINC 

Users are prompted for lower and upper limits of the vertical axis. 
YINC is the vertical axis increment between the ticmarks on the 
plot. 

Input 6a : PHMIN, PHMAX, PHINC (Twmplots only) 

For the program that generates twin plots (PSTWINPLOT.FOR), 
additional inputs for the vertical axis of the upper plot (usually the 
phase plot as in figure 2) are needed. PHMIN, PHMAX are the lower 
and upper limits of vertical axis, and PHINC is the coordinate 
increment between the ticmarks. 

Input 7 : TITLE 

Input title for plot, 60 characters maximum. However, to avoid 
overlapping problems, a maximum of 35 characters is recommended 
for the twinplot case. 
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Input 8 : Input first line of comments. 


Users can input two lines of comments. In addition, the program 
has an option to automatically put the date and time under the 
comment lines (see figures 1,2). This option is system dependent as 
can be seen on program line labeled 23 in subroutine PPLOT (program 
PSTWINPLOT.FOR) and in subroutine PPLOTXY (program 
PSPLOTXY.FOR). Thus the users are recommended to change the 
program line to their own system call. The program line is currently 
commented out to make the code portable. 

Input 3 : Input second line of comments 

Input 9a : Input next data file name (if multiple plots). 

Input 10 : PLOTLAB 

Input labels for the curves. Users are prompted to type in labels for 
each of the NPLOTS curves. 

The numerical labels at each of the ticmarks on the output plot is 
written with Fortran format F6.1. Program users can change this 
format to suit their need by changing appropriate parameters on 
program lines labeled 2002 and 2006 in subroutine PPLOT of 
PSPLOTXY.FOR and lines labeled 2002, 2006 and 2007 in subroutine 
PPLOT of PSTWINPLOT.FOR. Line labeled 2002 is for the horizontal 
axis. Line labeled 2006 is for formating numerical labels at the 
ticmarks of the vertical axis on the X-Y plot and for the vertical 
axis of amplitude plot on a twinplot. Line labeled 2007 is for the 
vertical axis of the phase plot in the twinplot case. 
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MAGNITUDE VS. ANGLE 

COMMENT UNE 1 
COMMENT UNE 2 
28-NOV-1990 08 49 30 47 



PLOT NUMBER 1 
PLOT NUMBER 2 
PLOT NUMBER 3 
PLOT NUMBER 4 



ANGLE IN DEGREES 


Figure 1: A sample X-Y plot (not true size) 



MAGNITUDE VS. ANGLE 

comment line 1 
comment line 2 
20-NOV-1 990 08 59 30 03 


Plot number 1 
Plot number 2 
Plot number 3 
Plot number 4 



ANGLE IN DEGREES 



Figure 2: A sample twinplot (not true size). 
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6.0 PROGRAM OUTPUT 


The programs generate postscript output files named either 
"RECPLOT.ps" for a X-Y plot program or "TWINPLOT.ps" for a twmplot 
program. 
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£★ ★ ★ ★ * ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ * ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ * 
C THIS PROGRAM GENERATES PLOTS IN XY COORDINATE. 

C THE OUTPUT FILE IS IN POSTSCRIPT LANGUAGE. 

C 

C -NUMERICAL LABELS AT THE TIC MARKS ARE WRITTEN 
C WITH FORMAT "F6.1" ON LINE LABEL 2002 AND 2006 IN 
C SUBROUTINE PPLOTXY. USER MAY WISH TO MODIFY THAT FORMAT 
C IF NEEDED . 

C 

C -MAXIMUM ARRAY SIZE IS 2000. TO INCREASE, CHANGE BOTH 
C PARAMETERS "NDIM" IN MAIN AND "NMAX" IN PPLOTXY 
C 

C -MAXIMUM 4 CURVES ON THE PLOT 
C 

C -LINE LABEL 23 IN SUBROUTINE PPLOTXY MAKES CALL TO A SYSTEM 
C DEPENDENT FUNCTION TO GET DATE AND TIME. USER SHOULD REPLACE 
C THIS LINE WITH A FUNCTION THAT WORKS ON THEIR SYSTEM. 

C FOR PORTABILITY THIS STATEMENT IS TEMPORARYLY COMMENTED OUT. 

C 1 

C -THE CONTENTS OF DATA FILES SHOULD BE IN ONE OF THE TWO 
C FOLLOWING FORMATS: 

C 1. FORMAT 1: 

C - First line : Npts, X(l), Xmc 

C where Npts = Number of Points, 

C X(l) = First data X value 

C Xinc = X increment (equally spaced data) 

C - Following lines: 

C Y (1) , PH (1) 

C Y (2) , PH (2) 

C Y (3) , PH (3) 

C ...... i 

C Where PH drray is a dummy array ('could be set to zero) 

C to maintain data files compability between PSPLOTXY.FOR 

C and PSPLOTAMPH . FOR 

C 

C 2. FORMAT 2: 

C - First line : Npts 

C where Npts = Number of Points, 

C - Following lines: 

C X(l), Y (1) , PH (1) 

C X(l), Y (2) , PH (2) 

C X(l), Y (3) , PH(3) 

C 

C Where PH array is a dummy array (could be set to zero) 

C to maintain data files compability between PSPLOTXY.FOR 

C and PSPLOTAMPH. FOR 

C 
C 

C 

PROGRAM PSPLOTXY 

C NPLOTS : NUMBER OF CURVES ON ONE PLOT 

PARAMETER (NDIM=2000) 

REAL* 4 AM (NDIM) , PH (NDIM) ,X (NDIM) 

WRITE (6, *) 

WRITE (6, 91) 'HOW MANY CURVES ON PLOT (4 MAX)? ' 
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READ (5, *) NPLOTS 

IF (NPLOTS . LE . 0 ) WRITE ( 6 ,*) ' ERROR NPLOTS ' ' 

C TYPES OF PLOTS: 

WRITE (6, *) 

WRITE (6,*) 'AXIS LABELS: ' 

WRITE (6, *) ' lbEF=0 : FREQ. SCAN (MAG VS FREQ) ' 

WRITE (6,*)' IDEF=1 : AZIM. SCAN (MAG VS ANGLE) ' 

WRITE (6,*)' IDEF=2 : PROMPT FOR AXIS LABELS ' 

WRITE (6,*)' IDEF=9 : PROMPT FOR PLOT SIZE & AXIS LABELS ' 

WRITE (6, *) ' (DEFAULT SIZE = 8.0 X 5.5 INCHES) ' 

WRITE (6, 91) 'Input IDEF: ' 

READ (5,*) IDEF 
WRITE (6, *) 

WRITE (6, 99) 

WRITE (6, 100) 

WRITE (6, 101) 

99 FORMAT (IX, 'TYPES OF INPUT DATA FORMAT (PH is dummy var) : ') 

100 FORMAT (IX, ' IFI = 1 :: header: N,XMIN,DX ; data: Y,PH') 

101 FORMAT (IX, ' =2 :: header: N ; data: X,Y,PH') 

91 FORMAT ( IX, A, $ ) 

WRITE (6, 91) 'Input IFI: ' 

READ (5,*) IFI 

DO 10 N=l, NPLOTS 

CALL DATAIN (NPTS, NDIM, X, AM, PH, IFI) 

CALL PPLOTXY (NPTS, NDIM, X, AM, PH, IDEF, NPLOTS, N, 3) 

10 CONTINUE 

STOP 

END 


CCC THIS IS INPUT READ FILE 
CCC 

SUBROUTINE DATAIN (NPTS, NDIM, X, AM, PH, IFI ) 
COMMON/HEADING/INFILE , LINE1 , LINE2, LINE3 
REAL X (NDIM) , AM (NDIM), PH (NDIM) 

CHARACTER* 60 LINE1, LINE2, LINE 3 
CHARACTER* 3 6 INFILE 
7 FORMAT (A36) 

91 FORMAT ( IX, A, $ ) 

WRITE (6, 91) 'ENTER (next) FILE NAME TO BE PLOTTED>. . . ' 
READ (5, 7) INFILE 

OPEN (UNIT=1 , NAME=INFILE, STATUS= ' OLD ' ) 

C FILE FORMAT READ: 

IF (IFI.EQ.l) THEN 

READ (1,*) NPTS, P_LO,P_INC i 

READ(1,*)( AM(I),PH(I), 1=1, NPTS ' ) 

DO 10 1=1, NPTS 

X (I) =P_LO+ (1-1 .0) *P_INC 
10 CONTINUE 

ELSE IF (IFI.EQ.2) THEN 
READ (1,*) NPTS 

READ (1, *) ( X(I),AM(I),PH(I), 1=1, NPTS ) 

ENDIF 

CLOSE (UNIT=1 , STATUS= ' KEEP ' ) 
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ooo o o O O 00000000000000 


RETURN 

END 


c 


SUBROUTINE PPLOTXY (NPTS, NDIM, X, AM, PH, IDEF, NPLOTS, NCOUNT, I GRID ) 


NPTS : INPUT NO. OF POINTS IN X AND Y ARRAY 
XX : INPUT REAL ARRAY X COORDINATES [1. .NPTS] 

AM : INPUT REAL ARRAY OF AMPLITUDES [1 . .NPTS] 

PH j INPUT REAL ARRAY OF PHASE DATA [1 . .NPTS] 

IDEF I INPUT INTEGER, IDEF=0 : FREQ. SCAN 

IDEF=1 : AZIM. SCAN 
IDEF=2 : CHANGE AXIS LABELS 
IDEF=9 : CHANGE X AND Y SIZE OF PLOT. 
TOTAL NUMBER OF CURVES ON PLOTS 
0: NO GRID (NOT INCORPORATED) 

1: HORIZONTAL GRID ONLY (NOT INCORPORATED) 

2: VERTICAL GRID ONLY (NOT INCORPORATED) 

3: BOTH HORIZONTAL AND VERTICAL GRID 
N th CURVE ON THE PLOT 1 


NPLOTS : 
IGRID : 


NCOUNT: 


COMMON/HEADING/ INFILE, LINE1, LINE2, LINE3 
PARAMETER (NMAX=2000) 

REAL X (NDIM) , AM (NDIM), PH (NDIM) 

REAL XX(NMAX) ,XAM(NMAX) ,XPH(NMAX) 

CHARACTER* 30 XLABEL, AMLABEL, PHLABEL 
CHARACTER* 30 ADASH ( 4 ) , PLOTLAB ( 4 ) , ATHICK 
CHARACTER* 60 LINE1, LINE2, LINE3, TITLE 
CHARACTER* 2 3 DATETIME 
CHARACTER* 3 6 INFILE 

SAVE XMIN, XMAX, YMIN, YMAX, AMMAX, AMMIN, PHMAX, PHMIN 
SAVE XLENGTH, YLENGTH 

CALL TO GET SYSTEM DATE & TIME: 

I IF (NCOUNT. EQ.l) CALL LIBSDATE TIME (DATETIME) 


SET DASH TYPES: 

ADASH (1)= 1 [) 0 setdash ’ 

ADASH (2 )= • [4 2] 0 setdash ' 

ADASH (3) = ' [6424] 0 setdash ' 

ADASH (4) = • [2 4] 0 setdash ' 

ATHICK® ' .6 setlinewidth ' 

TRANSFER THE VARIABLES, SO THE ROUTINE RETURNS THE ORIGINAL VARS 
DO 17 1=1, NPTS 
XX(I)=X(I) 

XAM(I) =AM(I) 

XPH(I)=PH(I) 

CONTINUE 


FOR MULTIPLE PLOT, GO TO 300 
IF (NCOUNT . GT . 1 ) GOTO 300 

SET UP DEFAULT PLOT SIZE IN INCHES: 
XLENGTH=5 .688 
YLENGTH=3 .438 
XLENGTH=8 . 000 
YLENGTH=5 . 500 
C SET UP AXIS LABELS: 
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IF (IDEF.EQ. 0) THEN 

XLABEL = ' FREQUENCY IN GHZ' 

AMLABEL= ' MAGNITUDE IN DB' 

ELSEIF ( IDEF . EQ . 1 ) THEN 

XLABEL = ' ANGLE IN DEGREES' 

AMLABEL= ' MAGNITUDE IN DB' 

ELSEIF (IDEF.EQ. 2) THEN 
WRITE (6, *) 

WRITE (6, 102) ' Enter the X label: ' 

READ (5, 101) XLABEL 

WRITE (6, 102) ' Enter the Y label: ' 

READ (5, 101) AMLABEL 
ELSEIF (IDEF.EQ. 9) THEN 
WRITE (6,*) 

WRITE (6, 102) ' Enter the physical size of X coord, (inch)' 
READ (5,*) XLENGTH 

WRITE (6, 102) ' Enter the physical size of Y coord, (inch)' 

READ (5,*) YLENGTH 

WRITE (6, 102) ' Enter the X label: ' 

READ (5, 101) XLABEL 

WRITE (6, 102) ' Enter the Y label: ' 

READ (5, 101) AMLABEL 
ELSE 

XLABEL = ' X ' 

AMLABEL= ' Y ' 

END IF 

101 FORMAT (30A) 

102 FORMAT ( IX, A, $ ) 

CX IF (XLENGTH. GT. 6) XLENGTH=6 . 0 

CX IF ( YLENGTH. GT. 4) YLENGTH=4 . 0 

IF (XLENGTH. GT. 8. 2) XLENGTH=8 . 2 
IF ( YLENGTH. GT. 6.0) YLENGTH=6. 0 
CALL AMINMAX (XX,NPTS,NDIM, AMAX,AMIN) 

WRITE (6,*) 

WRITE (6,*) ' XMIN= ', AMIN,' XMAX= ' , AMAX 

WRITE (6,*)' Enter Xmm, Xmax, Xinc for horizontal scaling ' 

READ (5,*) XMIN, XMAX, XINC 

NXTIC= INT ( (XMAX -XMIN) /XINC) 

CALL AMINMAX (XAM, NPTS, NDIM, AMAX, AMIN) 

WRITE (6,*) ' AMP_MIN= ', AMIN,' AMP_MAX= ' , AMAX 
WRITE(6,*) ' Enter Ymin, Ymax, Yinc for magnitude.' 

READ (5,*) AMMIN, AMMAX, AMINC 
NAMTIC= INT ( (AMMAX -AMMIN) /AMINC) 

100 FORMAT (A20) 

C END IF 

C CONTINUE HERE FOR MULTIPLE PLOTS 
300 CONTINUE 

C NORMALIZE THE ARRAY VALUES TO PHYSICAL DEVICE COORDS. 
TEMPX=XMAX-XMIN 
TEMP AM=AMMAX -AMMIN 
DO 1 1=1, NPTS 

IF (XX (I) . GT.XMAX) XX(I)=XMAX 
IF (XX (I) .LT.XMIN) XX (I) =XMIN 
XX(I)=( (XX (I) -XMIN) /TEMPX) * XLENGTH 
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IF (XAM(I) .GT.AMMAX) XAM(I)=AMMAX 
IF (XAM(I) .LT.AMMIN) XAM(I)=AMMIN 
XAM(I) = ( (XAM(I) -AMMIN) /TEMP AM) * YLENGTH 
XX(I)=XX(I)*72.0 
XAM(I) =XAM(I) *72.0 
1 CONTINUE 

XPHYINC = (XINC/TEMPX) *XLENGTH 
AMPHYINC = (AMINC/TEMPAM) * YLENGTH 

IF (NCOUNT . GT . 1 ) GOTO 305 

C TRANSFER TO PHYSICAL DEVICE COORDS MAX X AXIS 6 INCH MAX Y 7 INCH 

OPEN ( UNIT = 22, FILE = ' RECPLOT.ps 1 , 

+ RECL = 128, STATUS = 'UNKNOWN' ) 


I 


WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ’) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 


* % ' PS-Adobe ' 

'/fcopies 1 def' 

•/inch {72 mul} def' 

' 8.50 inch 0 inch translate' 

' 90 rotate ' 

' 1 inch 0.75 inch translate gsave' 
'/box{exch dup O' 

' rlineto exch dup 0 exch rlineto' 
' neg exch neg 0 rlineto' 

' 0 exch rlineto closepath) def' 

' .5 setlinewidth' 


WRITE (22, 1000) NXTIC 

1000 FORMAT (IX, '/nxtic ',13, ' def') 

WRITE (22, 1001) NAMTIC 

1001 FORMAT (IX, '/namtic ',13, ' def') 

WRITE (22, 1002) XLENGTH 

1002 FORMAT (IX, '/xsize ',f7.4, ' inch def') 
WRITE (22, 1003) YLENGTH 

1003 FORMAT (IX, '/ysize ',f7.4, ' inch def ' ) 
WRITE (22, 1004) XPHYINC 

1004 FORMAT (IX, '/xinc ',f7.4, ' inch def') 
WRITE (22, 1005) AMPHYINC 

1005 FORMAT (IX, '/amine ',f7.4, ' inch def') 
WRITE (22, 1006) NPTS 

1006 FORMAT (IX, ' / npts ',14, ' def') 

WRITE (22, 1007) DELY 

1007 FORMAT (IX, '/dely ',f7.3, ' inch def') 


& 


WRITE (22, FMT=' (A) ’ ) 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 

' xinc add exch 
WRITE (22, FMT=* (A) ') 


'/ticsize 7 def' 

'/do_x_tic (0 ticsize neg rlineto ' 

' 0 ticsize rmoveto) def ' 

'/do_x_axis {013-1 roll' 

' {do_x_tic currentpoint ' 

' exch xinc add exch moveto} for) def' 
'/do_y_tic {ticsize neg 0 rlineto ' 

' ticsize 0 rmoveto) def ' 

' /do_ver_grid { 1 sub 113-1 roll' 

' {currentpoint exch' , 

moveto currentpoint ysize add lineto ' 
' currentpoint ysize sub moveto}'. 
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& ' for stroke} def 

WRITE (22, FMT=' (A) ' ) ' /do_am_grid { 1 sub 1 1 3 -1 roll' 

WRITE (22, FMT=' (A) ' ) ' { currentpoint 
& ' amine add moveto currentpoint exch xsize add exch lineto' 

WRITE (22, FMT=' (A) ' ) ' currentpoint exch xsize suh ', 
fi 'exch moveto }for stroke} def' 

WRITE (22, FMT=' (A) ' ) ' /do_am_axis {013-1 roll' 

WRITE (22,FMT=' (A) ' ) ' {do_y_tic currentpoint' 

WRITE (22,FMT=' (A) ' ) ' amine add moveto}for} def' 

WRITE (22,FMT=’ (A) ’ ) ’ /do_alabel {currentpoint 3 -1 roll dup 1 
WRITE (22,FMT=' (A) ' ) ' stringwidth pop -2 div 0 rmoveto show' 
WRITE (22, FMT=' (A) • ) ' moveto }def' 

WRITE (22, FMT=' (A) ' ) ' /Cor_point { 4 . 5 inch xsize 2 div sub ' 
WRITE (22, FMT=' (A) ' ) ' .5 inch} def' 

WRITE (22, FMT=' (A) ' ) ' /Cor_point2{Cor_point ysize add 1.125 ' 
WRITE (22, FMT=' (A) ' ) * inch add} def' 

WRITE (22,FMT=' (A) ' ) '/Corner {Cor_point moveto} def' 

WRITE (22,FMT=' (A) ' ) ' /Corner2 {Cor_point2 moveto} def 1 

WRITE (22, FMT=' (A) ' ) ' Corner xsize ysize box stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner nxtic do_x_axis stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner 0 ysize ticsize add rmoveto ' 
WRITE (22, FMT=' (A) ' ) ' nxtic do_x_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner namtic do_am_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner xsize ticsize add 0 rmoveto' 
WRITE (22,FMT=' (A) ' ) ' namtic do_am_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner 0 ticsize -2.5 mul rmoveto ' 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 10 scalefont ' 

WRITE (22, FMT=' (A) ' ) 'setfont ' 

XTEMP=XMIN 
DO 1=1, NXTIC+1 
WRITE (22, 2002JXTEMP 

2002 FORMAT ( ' ( ' , F6. 1, ' ) do_alabel xinc 0 rmoveto ') 

XTEMP=XTEMP+XINC 

END DO 

WRITE (22,FMT=' (A) ' ) 'gsave Cor_point translate 90 rotate' 

WRITE (22, FMT=' (A) ' ) '0 0 moveto 0 ticsize 1.5 mul rmoveto ' 
XTEMP=AMMIN 
DO 1=1, NAMTIC+1 
WRITE (22, 2006)XTEMP 

2006 FORMAT ('(', F6 . 1, ' ) do_alabel amine 0 rmoveto ') 
XTEMP=XTEMP+AMINC 
END DO 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 16 scalefont ' 

WRITE (22, FMT=' (A) ' ) 'setfont 0 0 moveto' 

WRITE (22,FMT=' (A) ' ) 'ysize 2.0 div 3.5 ticsize mul rmoveto' 
WRITE (22, 2003) AMLABEL 

WRITE (22, FMT= ' (A) ' ) ' grestore ' 

WRITE (22,FMT=' (A) ' ) '/Helvetica findfont 16 scalefont setfont' 
WRITE (22,FMT=' (A) ' ) ' Corner xsize 2.0 div -5.5 ticsize mul ' 
WRITE (22, FMT= ' (A) ' ) 'rmoveto ' 

WRITE (22, 2003) XLABEL 

2003 FORMAT ( ' ( ' , A30 , ' ) do_alabel ' ) 
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305 


CONTINUE 


C PLOT LABEL, THREE COMMENT LINES & LINE STYLE MARKERS 
IF ( (NPLOTS • EQ . 1 ) .OR. (NCOUNT.EQ.NPLOTS) ) THEN 
WRITE (6, 91)' Enter TITLE (60 chars max): ' 

READ (5, 9) TITLE 
9 FORMAT (A60) 

91 FORMAT ( IX, A, $ ) 

WRITE (6,*) 'TWO LINES OF COMMENTS (DATE-TIME FOR LINE 3) : ' 

WRITE (6, 91) 'INPUT COMMENTS LINE 1>. . . ' 

READ (5, 9) LINE1 

WRITE (6, 91) ' COMMENTS LINE 2>. .. ' 

READ (5, 9) LINE2 

WRITE (6, 91)' COMMENTS LINE 3>. . .DATE-TIME ' 

C READ (5, 9) LINE 3 

WRITE (6, *) 

WRITE (22,FMT=' (A) ' ) ’/Helvetica findfont 16 scalefont setfont' 
WRITE (22, FMT=' (A) ') 'Corner -.6 inch ysize .8 inch add rmoveto' 

. WRITE (22, 2005) TITLE 

WRITE (22, FMT=' (A) ') '/Helvetica findfont 11 scalefont setfont' 
WRITE (22, FMT=' (A) ') 'Corner -.6 inch ysize .6 inch add rmoveto' 
WRITE (22, 2004) LINE1 

WRITE (22,FMT=' (A) ') 'Corner -.6 inch ysize .45 inch add rmoveto' 
WRITE (22, 2004) LINE2 

WRITE (22,FMT=' (A) ') 'Corner -.6 inch ysize .30 inch add rmoveto' 
WRITE (22, 2011) DATETIME 


C WRITE (22, 2004) LINE3 


2004 

FORMAT ( ' 

(', 

A60 

, ' ) show ' ) 

2005 

FORMAT (‘ 

(', 

A3 8 

, ' ) show ' ) 

2011 

FORMAT ( ' 

(', 

A23 

, ' ) show ' ) 


******** line STYLE MARKER ********* 


DO 111 1=1, NPLOTS 
WRITE (6, 113) I 

113 FORMAT (IX, 'Input LABEL for PLOT ',1X,I3, ',$) 

READ ( 5 , 1 1 4 ) PLOTLAB ( I ) 

114 FORMAT (A30) 

111 CONTINUE 

C WRITE THE PLOT LABEL 

WRITE (22, FMT=' (A) ') '/Helvetica findfont 11 scalefont setfont' 
DO 115 1=1, NPLOTS 
WRITE (22, 1201) ATHICK 
WRITE (22, 1201) AD ASH (I) 

IF (I.EQ.l)THEN 

WRITE (22,FMT=' (A) ') 'Corner 5.5 inch ysize .77 inch add rmoveto' 
ELSE IF (I. EQ. 2) THEN 

WRITE (22, FMT=' (A) ') 'Corner 5.5 inch ysize .62 inch add rmoveto' 
ELSE IF (I.EQ. 3) THEN 

WRITE (22, FMT=' (A) ') 'Corner 5.5 inch ysize .47 inch add rmoveto' 
ELSE IF (I.EQ. 4) THEN 

WRITE (22,FMT=' (A) ') 'Corner 5.5 inch ysize .32 inch add rmoveto' 
ENDIF 

WRITE (22, FMT=' (A) ') ' 0.5 inch 0 rlineto stroke' 

115 CONTINUE 

WRITE (22,FMT=' (A) ' ) '/Helvetica findfont 11 scalefont setfont' 
DO 117 1=1, NPLOTS 
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c 


WRITE (22, 1201) ADASH(I) 

IF ( I . EQ . 1 ) THEN 

WRITE (22,FMT=' (A) ') 'Corner 5.5 inch ysize .75 inch add rmoveto' 
ELSE IF (I. EQ. 2) THEN 

WRITE (22,FMT=' (A) ') 'Corner 5.5 inch ysize .60 inch add rmoveto' 
ELSE IF ( I . EQ . 3 ) THEN 

WRITE (22 / FMT=' (A) ') 'Corner 5.5 inch ysize .45 inch add rmoveto' 
ELSE IF (I. EQ. 4) THEN 

WRITE (22,FMT=' (A) ') 'Corner 5.5 inch ysize .30 inch add rmoveto' 
END IF 

WRITE (22, FMT=' (A) ' ) ' 0.4 inch 0 rmoveto' 

WRITE (22, 2005) PLOTLAB (I) 

117 CONTINUE 

c********end line style marker* ********* 

ENDIF 

WRITE (22, FMT=' (A) ' ) ' gsave Cor_point translate ' 

C TO SELECT DASH TYPES: 

WRITE (22, 1201) ATHICK 
WRITE (22, 1201) ADASH (NCOUNT) 

1201 FORMAT (IX, A30) 

WRITE (22, 3000) XX(1),XAM(1) 

3000 FORMAT (2 (F8 . 3, 2x) ' moveto ') 

WRITE (22, *) (XX(I) ,XAM(I) , ' lineto ',I=2,NPTS) 

WRITE (22, FMT = '(A)') 'stroke grestore' 

IF ( (NPLOTS.EQ.l) .OR. (NCOUNT. EQ.NPLOTS) ) THEN 
WRITE (22, FMT=' (A) ') '[1 5] 6 setdash .2 setlinewidth ' 

WRITE (22,FMT=' (A) ' ) ' Corner nxtic do_ver_grid ' 

WRITE (22,FMT=' (A) ' ) ' Corner namtic do_am_grid ' 

WRITE (22, FMT = '(A)') 'showpage' 

CLOSE ( 22 ) 

WRITE (6, *) 

WRITE (6,*) 'YOUR OUTPUT FILE IS "RECPLOT.ps" ' 

END IF 
RETURN 
END 

Finds min and max in an array of size NPTS 

SUBROUTINE AMINMAX (X, NPTS, NDIM, AMAX, AMIN) 

REAL X (NDIM) 

AMAX=-1000. 0 
AMIN=1000 . 0 

DO 1 1=1, NPTS , 

IF(X(I) .GT.AMAX) AMAX=X(I) 

IF (X (I) .LT.AMIN) AMIN=X ( I ) 

1 CONTINUE 

RETURN 
END 


32 



APPENDIX B 

Program PSTWINPLOT.FOR 
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Q* ********************************************************* ******** 

C THIS PROGRAM GENERATES TWO PLOTS ON THE SAME PAGE (AM VS X AND 
C PH VS X) . THE OUTPUT FILE IS IN POSTSCRIPT LANGUAGE . 

C 

C -NUMERICAL LABELS AT THE TIC MARKS ARE WRITTEN 
C WITH FORMAT "F6.1" ON LINE LABEL 2002, 2006 AND 2007 IN 
C SUBROUTINE PPLOT. USER MAY WISH TO MODIFY THAT FORMAT 
C IF NEEDED . 

C 

C -MAXIMUM ARRAY SIZE IS 1024. TO INCREASE, CHANGE BOTH 
C PARAMETERS "NDIM" IN MAIN AND "NMAX" IN PPLOT 
C 

C -MAXIMUM 4 CURVES ON THE PLOT 
C 

C -LINE LABEL 23 IN SUBROUTINE PPLOT MAKES CALL TO A SYSTEM 
C DEPENDENT FUNCTION TO GET DATE AND TIME. USER SHOULD REPLACE 
C THIS LINE WITH A FUNCTION THAT WORKS ON THEIR SYSTEM. 

C FOR. PORTABILITY THIS STATEMENT IS TEMPORARYLY COMMENTED OUT. 

C 

C -THE CONTENTS OF DATA FILES SHOULD BE IN ONE OF THE TWO 
C FOLLOWING FORMATS : 

C 1. FORMAT 1: 

C - First line : Npts, X(l), Xinc 

C where Npts = Number of Points, 

C X (1) = First data X value 

C Xinc = X increment (equally spaced data) 

C - Following lines: 

C Y (1) , PH (1) 

C Y (2) , PH (2) 

C Y (3) , PH (3) 

C 

c 

C 2. FORMAT 2: 

C - First line : Npts 

C where Npts = Number of Points, 

C - Following lines: 

C X (1) , Y (1) , PH (1) 

C X(l), Y (2) , PH (2) 

C X(l), Y (3) , PH (3) 

C 

c 

c 

PROGRAM PSPLOT 

Plot Amp. and Phase data to the postscript printer 
NPLOTS : NUMBER OF CURVES ON ONE PLOT 
PARAMETER (NDIM=1024) 

REAL* 4 AM (NDIM) , PH (NDIM) ,X (NDIM) 

WRITE (6, *) 

WRITE (6,*) 'THIS PROGRAM PLOTS TWIN PLOTS (MAG & PHASE)' 
WRITE (6, 91) 'HOW MANY CURVES ON PLOT (4 MAX)? ' 

READ (5,*) NPLOTS 

IF (NPLOTS. LE.0) WRITE ( 6 ,*)' ERROR NPLOTS ! ' 

C TYPES OF PLOTS: 

WRITE (6, *) 

WRITE (6,*) 'AXIS LABELS: ' 
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WRITE (6,*) ' 

IDEF=0 : 

WRITE (6, * ) ' 

IDEF=1 : 

WRITE (6,*) ' 

IDEF=2 : 

WRITE (6,*) ' 

IDEF=9 : 

WRITE (6, *) ' 



XLENGTH=5 . 688 
YLENGTH=3 .438 


FREQ. SCAN (MAG VS FREQ) ' 

AZIM. SCAN (MAG VS ANGLE) ' 

PROMPT FOR AXIS LABELS ' 

PROMPT FOR PLOT SIZE & AXIS LABELS ' 
(DEFAULT SIZE= 5.688 X 3.438 INCHES)' 


WRITE (6, 91) 'Input IDEF: ' 

READ (5,*) IDEF 
WRITE (6,*) 

WRITE (6, 99) 

WRITE (6, 100) 

WRITE (6, 101) 

99 FORMAT (IX, 'TYPES OF INPUT DATA FORMAT : ') 

100 FORMAT (IX, * IFI = 1 :: header: N,XMIN,DX ; data: Y,PH') 

101 FORMAT (IX, ' = 2 :: header: N ; data: X,Y,PH') 

91 FORMAT ( IX, A, $ ) 

WRITE (6, 91) 'Input IFI: ' 

READ (5, *) IFI 

DO 10 N=l, NPLOTS 

CALL DATAIN (NPTS, NDIM, X, AM, PH, IFI) 

CALL PPLOT (NPTS, NDIM, X, AM, PH, IDEF, NPLOTS, N, 3) 

10 CONTINUE 


STOP 

END 


CCC THIS IS INPUT READ FILE 
CCC 

SUBROUTINE DATAIN (NPTS, NDIM,X, AM, PH, IFI) 

COMMON/ HEADING/ INFILE , LINE1, LINE2, LINE 3 
PARAMETER (NMAX=1024) 

REAL X (NDIM) , AM (NDIM) , PH (NDIM) 

CHARACTER* 60 LINE1, LINE2, LINE3 
CHARACTER* 36 INFILE 
7 FORMAT (A36) 

91 FORMAT ( IX, A, $ ) 

WRITE (6, 91) 'ENTER (next) FILE NAME TO BE PLOTTED>. . . ' 
READ (5, 7) INFILE 

OPEN (UNIT=1 , NAME=INFILE, STATUS= ' OLD ' ) 


C FILE FORMAT READ: 

IF (IFI.EQ.l) THEN 

READ (1, *) NPTS, P_LO, P_INC 

READ (1, *) ( AM(I) , PH (I) , 1=1, NPTS ) 

DO 10 1=1, NPTS 

X ( I ) =P_LO+ ( 1-1 . 0 ) *P_INC 
10 CONTINUE 

ELSE IF (IFI.EQ.2) THEN 
READ (1,*) NPTS 

READ(1,*)( X(I),AM(I),PH(I), 1=1, NPTS ) 
ENDIF 

CLOSE (UNIT=1 , STATUS= ' KEEP ' ) 


RETURN 

END 
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c 

SUBROUTINE PPLOT (NPTS, NDIM, X, AM, PH, IDEF, NPLOTS, NCOUNT, I GRID) 

NPTS : INPUT NO. OF POINTS IN X AND Y ARRAY 

XX : INPUT REAL ARRAY X COORDINATES [1 . .NPTS] 

AM : INPUT REAL ARRAY OF AMPLITUDES [1 . .NPTS] 

PH : INPUT REAL ARRAY OF PHASE DATA [ 1 . . NPTS ) 

IDEF : INPUT INTEGER, IDEF=0 : FREQ. SCAN 

IDEF=1 : AZIM. SCAN 
IDEF=2 : CHANGE AXIS LABELS 
IDEF=9 : CHANGE X AND Y SIZE OF PLOT. 
NPLOTS : TOTAL NUMBER OF CURVES ON PLOTS 

IGRID : 0: NO GRID (NOT INCORPORATED) 

1: HORIZONTAL GRID ONLY (NOT INCORPORATED) 

2: VERTICAL GRID ONLY (NOT INCORPORATED) 

3: BOTH HORIZONTAL AND VERTICAL GRID 
NCOUNT: N th CURVE ON THE PLOT 

COMMON/HEADING/INFILE, LINE1, LINE2, LINE3 
PARAMETER (NMAX=1024) 

REAL X (NDIM) , AM (NDIM), PH (NDIM) 

REAL XX(NMAX) ,XAM(NMAX) ,XPH(NMAX) 

CHARACTER* 30 XLABEL, AMLABEL, PHLABEL 
CHARACTER* 3 0 ADASH (4) , PLOTLAB (4 ) 

CHARACTER* 60 LINE1, LINE2, LINE3, TITLE 
CHARACTER* 2 3 DATETIME 
CHARACTER* 3 6 INFILE 

SAVE XMIN, XMAX, YMIN, YMAX, AMMAX, AMMIN, PHMAX, PHMIN 
SAVE XLENGTH, YLENGTH 

C CALL TO SYSTEM TO GET DATE & TIME 
23 IF (NCOUNT. EQ.l) CALL LIB$DATE_TIME (DATETIME) 

C SET DASH TYPES: 

ADASH (1)= ' [] 0 setdash ' 

ADASH (2) = ' [4 2] 0 setdash ' 

ADASH ( 3 ) = ' [6 4 2 4] 0 setdash ' 

ADASH (4) = ' [2 4] 0 setdash ' ; 

C TRANSFER THE VARIABLES, SO THE ROUTINE RETURNS THE ORIGINAL VARS 
DO 17 1=1, NPTS 
XX(I)=X(I) 

XAM ( I ) =AM ( I ) 

XPH (I) =PH (I) 

17 CONTINUE 

IF (NCOUNT. GT.l) GOTO 300 

CXX 

C SET UP DEFAULT PLOT SIZE IN INCHES: 

XLENGTH=5 .688 
YLENGTH=3. 438 
XLENGTH=8 .000 
YLENGTH=5 .500 
SET UP AXIS LABELS: 

IF (IDEF.EQ. 0) THEN 

XLABEL = • FREQUENCY IN GHZ' 

AMLABEL= ' MAGNITUDE IN DB' 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 
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PHLABEL= ' PHASE IN DEGREES' 

ELSEIF ( IDEF . EQ . I ) THEN 

XLABEL = ' ANGLE IN DEGREES' 

AMLABEL= ' MAGNITUDE IN DB' 

PHLABEL= ' PHASE IN DEGREES' 

ELSEIF (IDEF.EQ.2) THEN 
WRITE (6, *) 

WRITE (6, 102) ' Enter the X label: ' 

READ (5, 101) XLABEL 

WRITE (6, 102) ' Enter the AM (amplitude) label: ' 

READ (5, 101) AMLABEL 

WRITE (6, 102) ' Enter the PH (phase) label: ' 

READ (5, 101) PHLABEL 
ELSEIF (IDEF.EQ. 9) THEN 
WRITE (6, *) 

WRITE (6, 102) ' Enter the physical size of X coord, (inch)' 
READ (5,*) XLENGTH 

WRITE (6, 102) ' Enter the physical size of Y coord, (inch)' 

READ (5,*) YLENGTH 

WRITE (6, 102) ' Enter the X label: ' 

READ (5, 101) XLABEL 

WRITE (6, 102) ' Enter the AM (amplitude) label: ' 

READ (5, 101) AMLABEL 

WRITE (6, 102) ' Enter the PH (phase) label: ' 

READ (5, 101) PHLABEL 
ELSE 

XLABEL = ' X ' 

AMLABEL= ' MAGNITUDE IN DB ' 

PHLABEL= ' PHASE IN DEGREES' 

END IF 

101 FORMAT (30A) 

102 FORMAT ( IX, A, $ ) 

IF (XLENGTH. GT. 6) XLENGTH=6 . 0 
IF ( YLENGTH. GT. 4) YLENGTH=4 . 0 
CALL AMINMAX (XX, NPTS, NDIM, AMAX, AMIN) 

WRITE (6,*) 

WRITE (6,*) ' XMIN= ', AMIN,' XMAX= ' , AMAX 

WRITE (6,*) 'Enter Xmin, Xmax, Xinc for horizontal scaling' 

READ (5,*) XMIN, XMAX, XINC 

NXTIC= INT ( (XMAX-XMIN) /XINC) 

CALL AMINMAX (XAM, NPTS, NDIM, AMAX, AMIN) 

WRITE (6,*) ' AMP_MIN= ', AMIN,' AMP_MAX= ' , AMAX 

WRITE (6, *) ' Enter min, max, inc for magnitude.' 

READ (5, *) AMMIN, AMMAX, AMINC 
NAMTIC= INT ( (AMMAX -AMMIN) /AMINC) 

CALL AMINMAX (XPH, NPTS, NDIM, AMAX, AMIN) 

WRITE (6,*) ' PH_MIN= ', AMIN,' PH_MAX= ' , AMAX 

WRITE (6, *) ' Enter min, max, inc for phase.' 

READ (5, *) PHMIN, PHMAX, PHINC 
NPHTIC= INT ( (PHMAX-PHMIN) /PHINC) 

100 FORMAT (A20) 

C CONTINUE HERE FOR MULTIPLE PLOTS 
300 CONTINUE 

C NORMALIZE THE ARRAY VALUES TO PHYSICAL DEVICE COORDS. 
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TEMPX=XMAX-XMIN 
TEMPAM=AMMAX-AMMIN 
TEMPPH=PHMAX-PHMIN 
DO 1 1=1, NPTS 

IF (XX (I) .GT.XMAX) XX(I)=XMAX 
IF (XX(I) .LT.XMIN) XX (I) =XMIN 
XX (I) = ( (XX (I) -XMIN) /TEMPX) * XLENGTH 
IF (XAM(I) . GT.AMMAX) XAM(I)=AMMAX 
IF (XAM(I) .LT.AMMIN) XAM(I)=AMMIN 
XAM (I) = ( (XAM(I) -AMMIN) /TEMP AM) * YLENGTH 
IF (XPH(I) .GT.PHMAX) XPH(I)=PHMAX 
IF (XPH(I) .LT.PHMIN) XPH(I)=PHMIN 
XPH (I) = ( (XPH (I) -PHMIN) /TEMPPH) * YLENGTH 
XX(I)=XX(I)*72.0 
XAM(I)=XAM(I)*72.0 
XPH(I)=XPH(I)*72.0 
1 CONTINUE 

XPHYINC = (XINC/TEMPX) * XLENGTH 
AMPHYINC = (AMINC/TEMPAM) * YLENGTH 
PHPHYINC = (PHINC/ TEMPPH) *YLENGTH 

IF (NCOUNT . GT . 1 ) GOTO 305 

C TRANSFER TO PHYSICAL DEVICE COORDS MAX X AXIS 6 INCH MAX Y 7 INCH 

OPEN ( UNIT = 22, FILE = 'TWINPLOT.ps', 

+ STATUS = 'UNKNOWN' ) 


WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE ( 22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ') 
WRITE (22, FMT=' (A) ' ) 
WRITE (22, FMT=' (A) ' ) 


' % ' PS-Adobe ' 

'/#copies,l def' 

'/inch {72 mul) ! def' 

• 1 inch 1 inch translate gsave' 
'/box(exch dup O' 

' rlineto exch dup 0 exch rlineto' 
' neg exch neg 0 rlineto' 

' 0 exch rlineto closepath) def' 

' .5 setlinewidth' 


WRITE (22, 1000) NXTIC 

1000 FORMAT (IX, '/nxtic ',13, ' def') 

WRITE (22, 1001) NAMTIC 

1001 FORMAT (IX, '/namtib ',13, ' def') 

WRITE (22, 1009) NPHTIC 

1009 FORMAT (IX, '/nphtic ',13, ' def') 

WRITE (22, 1002) XLENGTH 

1002 FORMAT (IX, '/xsize ',f7.4, ' inch def') 
WRITE (22, 1003) YLENGTH 

1003 FORMAT (IX, '/ysize ',f7.4, • inch def') 
WRITE (22, 1004) XPHYINC 

1004 FORMAT (IX, '/xinc ',f7.4, ' inch def') 
WRITE (22, 1005) AMPHYINC 

1005 FORMAT (IX, '/amine ',f7.4, ' inch def') 
WRITE (22, 1008) PHPHYINC 

1008 FORMAT (IX, '/phinc ',f7.4, ' inch def') 

WRITE (22, 1006) NPTS 

1006 FORMAT (IX, '/npts ',14, • def') 

WRITE (22, 1007) DELY 
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1007 FORMAT (IX, '/dely ', 11 . 2 , ' inchdef') 

WRITE ( 22 , FMT= ' (A) ' ) '/ticsize 7 def' 

WRITE (22, FMT=' (A) ' ) ' /do_x_tic {0 ticsize neg rlineto 1 
WRITE (22,FMT= f =' (A) ' ) 1 0 ticsize rmoveto} def ' 

WRITE (22, FMT=' (A) ') ' /do_x_axis (013-1 roll' 

WRITE (22, FMT=' (A) ' ) ' (do_x_tic currentpoint ' 

WRITE (22, FMT=' (A) ') ' exch xinc add exch moveto}for) def' 

WRITE (22,FMT=' (A) ' ) ' /do_y_tic (ticsize neg 0 rlineto ' 

WRITE (22, FMT=' (A) ') ' ticsize 0 rmoveto} def ' 

WRITE (22, FMT=' (A) ' ) ' /do_ver_grid { 1 sub 1 1 3 -1 roll' 

WRITE (22,FMT=' (A) ' ) ' (currentpoint exch' , 

& ' xinc add exch moveto currentpoint ysize add lineto ' 

WRITE (22,FMT=' (A) ' ) ' currentpoint ysize sub moveto}', 

& ' for stroke} def' 

WRITE (22, FMT=' (A) ' ) '/do_am_grid { 1 sub 1 1 3 -1 roll' 

WRITE (22, FMT=' (A) ' ) ' (currentpoint ' , 

& ' amine add moveto currentpoint exch xsize add exch lineto' 

WRITE (22,FMT=' (A) ' ) ' currentpoint exch xsize sub ', 

& 'exch moveto }for stroke} def' 

WRITE (22, FMT=' (A) ' ) • /do_ph_grid { 1 sub 1 1 3 -1 roll' 

WRITE (22, FMT=' (A) ' ) ' (currentpoint ' , 

£ ' phinc add moveto currentpoint exch xsize add exch lineto' 

WRITE (22,FMT=' (A) ' ) ' currentpoint exch xsize sub ', 

& 'exch moveto }for stroke} def' 

WRITE (22, FMT=' (A) ' ) ' /do_am_axis (013-1 roll' 

WRITE (22, FMT=' (A) ' ) ' (do_jr_tic currentpoint' 

WRITE (22, FMT=' (A) ' ) ' amine add moveto}for} def' 

WRITE (22, FMT=' (A) ') ' /do_ph_axis (013-1 roll' 

WRITE (22,FMT=' (A) ' ) ' (do_y_tic currentpoint' 

WRITE (22, FMT=' (A) ') ' phinc add moveto}for} def' 

WRITE (22,FMT=' (A) ' ) ' /do_alabel {currentpoint 3 -1 roll dup ' 
WRITE (22, FMT=' (A) ' ) ' stnngwidth pop -2 div 0 rmoveto show' 
WRITE (22,FMT=' (A) ' ) ' moveto }def' 

WRITE (22, FMT=' (A) ' ) ' /Cor_point{ 3. 5 inch xsize 2 div sub ' 
WRITE (22, FMT=' (A) ' ) ’ .5 inch} def' 

WRITE (22,FMT=' (A) ' ) ' /Cor_point2(Cor_point ysize add 1.125 ' 
WRITE (22, FMT=' (A) ') ' inch add} def' 

WRITE (22,FMT=' (A) ' ) '/Corner (Cor_point moveto} def' 

WRITE (22,FMT=' (A) ' ) ' /Corner2 (Cor_point2 moveto} def ' 

WRITE (22, FMT=' (A) ') ' Corner xsize ysize box stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner nxtic do_x_axis stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner 0 ysize ticsize add rmoveto ' 
WRITE (22, FMT=' (A) ') ' nxtic do_x_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner namtic do_am_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner xsize ticsize add 0 rmoveto' 
WRITE (22, FMT=' (A) ' ) ' namtic do_am_axis stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner2 xsize ysize box stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner2 nxtic do_x_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner2 0 ysize ticsize add rmoveto ' 
WRITE (22, FMT=' (A) ') ' nxtic do_x_axis stroke' 

WRITE (22, FMT=' (A) ' ) ' Corner2 nphtic do_ph_axis stroke' 

WRITE (22,FMT=' (A) ' ) ' Corner2 xsize ticsize add 0 rmoveto' 
WRITE (22,FMT=' (A) ' ) ' nphtic dojph_axis stroke' 
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WHITE(22, FMT=' (A) ' ) ' Corner 0 ticsize -2.5 mul rmoveto ' 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 10 scalefont ' 

WRITE (22, FMT=' (A) ' ) 'setfont ' 

XTEMP=XMIN 
DO 1=1, NXTIC+1 
WRITE (22,2002)XTEMP 

2002 FORMAT (' (',F6.1, ') do_alabel xinc 0 rmoveto ') 
XTEMP=XTEMP+XINC 

END DO 

WRITE (22, FMT=' (A) ' ) ' Corner2 0 ticsize -2.5 mul rmoveto ' 
XTEMP=XMIN 
DO 1=1, NXTIC+1 
WRITE (22, 2002JXTEMP 
XTEMP=XTEMP+XINC 
END DO 

WRITE (22,FMT=‘ (A) ' ) 'gsave Cor_point translate 90 rotate' 
WRITE (22,FMT=‘ (A) ' ) '0 0 moveto 0 ticsize 1.5 mul rmoveto ' 
XTEMP=AtWIN 
DO 1=1, NAMTIC+1 
WRITE (22,2006)XTEMP 

2006 FORMAT (' (' ,F6.1, ' ) do_alabel amine 0 rmoveto ') 
XTEMP=XTEMP+AMINC 

END DO 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 16 scalefont ' 

WRITE (22, FMT=' (A) ') 'setfont 0 0 moveto' 

WRITE (22,FMT=' (A) ' ) 'ysize 2.0 div 3.5 ticsize mul rmoveto' 
WRITE (22, 2003) AMLABEL 

WRITE (22, FMT=' (A) ' ) 'grestore Cor_point2 translate 90 rotate' 
WRITE (22, FMT=' (A) ' ) '0 0 moveto 0 ticsize 1.5 mul rmoveto ' 
XTEMP=PHMIN 

DO 1=1, NPHTIC+1 ' 

WRITE (22, 2007)XTEMP 

2007 FORMAT ( ' ( ' , F6 . 1, ' ) do_alabel phinc 0 rmoveto ') 
XTEMP=XTEMP+PHINC 

END DO 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 16 scalefont ' 

WRITE (22, FMT=' (A) ') 'setfont 0 0 moveto' 

WRITE (22, FMT=' (A) ' ) 'ysize 2.0 div 3.5 ticsize mul rmoveto' 

WRITE (22, 2003) PHLABEL 

WRITE (22, FMT=' (A) ' ) ' grestore ' 

WRITE (22, FMT=' (A) ' ) '/Helvetica findfont 16 scalefont setfont' 
WRITE (22, FMT=' (A) ' ) ' Corner xsize 2.0 div -5.5 ticsize mul ' 
WRITE (22, FMT= ' (A) ' ) ' rmoveto ' 

WRITE (22, 2003) XLABEL 

2003 FORMAT (• (', A30 ,' ) do_alabel ') 

WRITE (22, FMT=* (A) ' ) ' Corner2 xsize 2.0 div -5.5 ticsize mul ' 
WRITE (22, FMT=' (A) ' ) ' rmoveto ' 

WRITE (22, 2003) XLABEL 

305 CONTINUE 

C PLOT LABEL, THREE COMMENT LINES & LINE STYLE MARKERS 
IF( (NPLOTS . EQ . 1 ) .OR. (NCOUNT.EQ.NPLOTS) ) THEN 
WRITE (6,*) 

WRITE (6, 91)' Enter TITLE (40 chars max): ' 
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READ (5, 9) TITLE 
9 FORMAT (A60) 

91 FORMAT ( IX, A, $ ) 

WRITE (6,*) 'TWO LINES OF COMMENTS (DATE-TIME FOR LINE 3):' 

WRITE (6, 91) 'INPUT COMMENTS LINE 1>. . . ' 

READ (5, 9) LINE1 

WRITE (6, 91) ' COMMENTS LINE 2>. . . ' 

READ (5, 9) LINE? 

WRITE (6, 91)' COMMENTS LINE 3>. . .DATE-TIME ' 

C READ (5, 9) LINE 3 

WRITE (6, *) 

WRITE (22,FMT=' (A) ' ) '/Helvetica findfont 16 scalefont setfont" 
WRITE (22,FMT=' (A) ' ) 'Corner2 -.6 inch ysize .8 inch add rmoveto' 
WRITE (22, 2005) TITLE 

WRITE (22,FMT=' (A) ' ) '/Helvetica findfont 11 scalefont setfont' 
WRITE (22,FMT=' (A) ' ) 'Corner2 -.6 inch ysize .6 inch add rmoveto' 
WRITE (22, 2004) LINE1 

WRITE (22, FMT=' (A) ' ) 'Corner2 -.6 inch ysize .45 inch add rmoveto' 
WRITE ( 22 , 2 0 0 4 ) LINE2 

WRITE (22,FMT=' (A) ' ) 'Corner2 -.6 inch ysize .30 inch add rmoveto* 
C WRITE (22, 2004) LINE3 

WRITE (22, 2011) DATETIME 


2004 

FORMAT ( ' 

(\ 

A60 

, ' ) show ' 

') 

2005 

FORMAT ( ' 


A38 

, ' ) show ' 

') 

2011 

FORMAT ( ' 

<\ 

A23 

, ' ) show 

') 


C******** LINE STYLE MARKER ********* 

C 

DO 111 1=1, NPLOTS 
WRITE (6, 113) I 

113 FORMAT (IX, 'Input LABEL for PLOT ',1X,I3, * : ',$) 

READ (5, 114 ) PLOTLAB (I ) 

114 FORMAT (A30) 

111 CONTINUE 

C WRITE THE PLOT LABEL 

WRITE (22,FMT=' (A) ' ) '/Helvetica findfont 11 scalefont setfont' 
DO 115 1=1, NPLOTS 
WRITE (22, 1201) ADASH(I) 

IF (I . EQ. 1) THEN 

WRITE (22,FMT=' (A) ' ) 'Corner2 3.5 inch ysize .77 inch add rmoveto' 
ELSE IF (I. EQ. 2) THEN 

WRITE (22,FMT=' (A) ' ) 'Corner2 3.5 inch ysize .62 inch add rmoveto' 
ELSE IF (I. EQ. 3) THEN 

WRITE (22, FMT=' (A) ') 'Corner2 3.5 inch ysize .47 inch add rmoveto' 
ELSE IF (I. EQ. 4) THEN 

WRITE (22, FMT=' (A) ') 'Corner2 3.5 inch ysize .32 inch add rmoveto' 
ENDIF 

WRITE (22, FMT=' (A) ' ) ' 0.5 inch 0 rlineto stroke' 

115 CONTINUE 

WRITE (22, FMT= ' (A) ' ) '/Helvetica findfont 11 scalefont setfont' 
DO 117 1=1, NPLOTS 
WRITE (22, 1201) AD ASH (I) 

IF (I .EQ. 1) THEN 

WRITE (22, FMT=' (A) ') 'Corner 2 3.5 inch ysize .75 inch add rmoveto' 
ELSE IF (I.EQ.2)THEN 

WRITE (22,FMT=* (A) ' ) 'Corner2 3.5 inch ysize .60 inch add rmoveto' 
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ELSE IF (I.EQ. 3) THEN 

WRITE (22,FMT t =' (A) ') 'Corner 2 3.5 inch ysize .45 inch add rmoveto' 
ELSE IF (I.EQ. 4) THEN 

WRITE (22, FMT=' (A) ' ) 'Corner2 3.5 inch ysize .30 inch add rmoveto' 
ENDIF 

WRITE (22, FMT= ' (A) ' ) ' 0.4 inch 0 rmoveto' 

WRITE (22, 2005) PLOTLAB (I ) 

117 CONTINUE 


c********end line style marker* ********* 

ENDIF 

WRITE (22, FMT=' (A) ' ) ' gsave Cor_point translate ' 

C TO SELECT DASH TYPES: 

WRITE (22, 1201) ADASH (NCOUNT) 

1201 FORMAT (IX, A30) 

WRITE (22, 3000) XX(1),XAM(1) 

3000 FORMAT (2 (F8. 3, 2x) ' moveto ') 

WRITE (22, *) (XX(I) ,XAM(I) , ' lineto ',I=2,NPTS) 

WRITE (22, FMT = ' (A) ' ) ' stroke grestore ' 

WRITE (22, FMT= ' (A) ' ) ' gsave Cor_point2 translate ' 
C TO SELECT DASH TYPES: 

WRITE (22, 1201) ADASH (NCOUNT) 

WRITE (22, 3000) XX(1),XPH(1) 

WRITE (22, *) (XX(I) ,XPH(I) , ' lineto • ,I=2,NPTS) 
WRITE (22, FMT = '(A)') 'stroke grestore ' 


IF ( (NPLOTS . EQ . 1 ) . OR . (NCOUNT . EQ . NPLOTS ) ) THEN 

WRITE (22, FMT=' (A) ' ) '[1 5] 6 setdash .2 setlinewidth ' 

WRITE (22, FMT=' (A) ' ) ' Corner nxtic do_ver_grid ' 

WRITE (22,FMT=' (A) ' ) ' Corner2 nxtic do_ver_grid ' 

WRITE (22,FMT=' (A) ' ) ' Corner namtic do_am_grid ' 

WRITE (22, FMT=' (A) ') ' Corner2 nphtic do__ph_grid' 

WRITE (22, FMT = '(A)') 'showpage' 

CLOSE ( 22 ) 

WRITE (6,*) 

WRITE (6,*) 'THE OUTPUT FILE IS "TWINPLOT.ps" ' 

END IF i 

RETURN 
END 

C 

C Finds min and max in an array of size NPTS 

SUBROUTINE AMINMAX (X, NPTS, NDIM,AMAX, AMIN) 

REAL X (NDIM) 

AMAX=-1000.0 
AMIN=1000 . 0 
DO 1 1=1, NPTS 

IF(X(I) .GT.AMAX) AMAX=X(I) 

IF(X(I) .LT.AMIN) AMIN=X ( I ) 

1 CONTINUE 

RETURN 
END 
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